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O ■ ABSTRACT 
(N 

^ I Two-dimensional particle simulations of a relativistic Harris current sheet of 

pair (e^) plasmas have demonstrated that the system is unstable to the rela- 
'nI" I tivistic drift kink instability and that a new kind of acceleration process takes 

place in the deformed current sheet. This process contributes to the generation 
^ I of non-thermal particles and contributes to the fast magnetic dissipation in the 

I current sheet structure. The acceleration mechanism and a brief comparison with 



relativistic magnetic reconnection are presented. 



I Subject headings: acceleration of particles — instabilities; plasmas — relativity; 



magnetic fields — pulsars (individual: Club Pulsar) 
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^ , 1. Introduction 

A current sheet structure with reversed magnetic fields can be seen everywhere in the 
^ universe. Among the magnetic field dissipation processes in the current sheet, magnetic 

' reconnection may be the most popular one. Breaking the topology of magnetic field lines, 
it releases the magnetic energy into particles' kinetic energy in a relatively short time-scale 
(several to several tens of an Alfven transit time ta = A/Va, where A is the half thickness 
of the current sheet and Va is the typical Alfven speed). It plays an essential role in solar 
terrestrial sites such as the Earth's magnetotail (Dungey 1961) and source regions of the solar 
flare (Parker 1957). Moreover, it is expected in astrophysical sites such as active galactic 
nuclei (Lesch & Birk 1998; Schopper et al. 1998) and pulsar winds, especially their well- 
studied example of the Crab Nebula (Coroniti 1990). In case of the Crab Nebula, energy 
conversion in a relativistic wind of plasmas has been a long-standing problem (the so- 
called the "cT-problem" , where a is the ratio of the Poynting flux energy to the particle kinetic 
flux) (Kennel & Coroniti 1984; Coroniti 1990; Lyubarsky & Kirk 2001; Kirk & Skj^raasen 
2003). The wind is originally Poynting-dominated (a ~ 10^) close to the neutron star 
(Arons 1979); but it is kinetic-dominated near the termination shock (a is far less than 
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the unity; a ~ 10"^ by Kennel & Coroniti (1984) and a ~ 10"^ by Mori et al. (2004)). 
Recently, relativistic magnetic reconnection in the striped current sheets has been discussed 
as the most possible dissipation mechanism (Coroniti 1990; Lyubarsky & Kirk 2001; Kirk 
&; Skjasraasen 2003). The properties of relativistic reconnection in plasmas have been 
investigated by an analytical study (Blackman & Field 1994) and by computer simulations 
(Zenitani & Hoshino 2001; Jaroschek et al. 2004). One of the most important features of 
relativistic reconnection is the generation of non-thermal particles due to direct acceleration 
by dc electric field around the X-type region (Zenitani & Hoshino 2001). However, the above 
discussion is based on a two-dimensional picture of reconnection. Instabilities in a cross-field 
plane are also of critical importance to understanding a realistic current sheet problem 
that includes reconnection. Does reconnection grow faster than cross-field instabilities? Is 
reconnection the fastest dissipation process? In the case of non-relativistic ion-electron 
plasmas, the following instabilities have been discussed: the lower hybrid drift instability 
(LHDl) (Krall 1971; Davidson et al. 1977), the Kelvin- Helmhohz instability (KHl) (Yoon & 
Drake 1996; Shinohara et al. 2001), and the drift kink instability (DKl) (Zhu & Winglee 1996; 
Pritchett et al. 1996; Daughton 1998). They lead to plasma heating and particle acceleration 
as well as the triggering of reconnection. The purpose of this Letter is to search for the basic 
properties of such cross-field instabilities and of the non-thermal particle generation in a 
relativistic current sheet. In this Letter, we report our 2D simulations in a cross-field 
plane of a relativistic current sheet and we report a new acceleration process related to 
the relativistic drift kink instability (RDKI) (Zenitani & Hoshino 2003). 



2. Simulation 



The simulation code used here is a three-dimensional particle-in-cell code. The system 
size is composed of 1 (X) x 256 (Y) x 512 [Z) grids. We consider periodic boundaries in 
the X, Y directions, and we periodically set two simulation boxes in the Z direction. The 
typical scale of the current sheet A is set to 10 grids, so that the Y boundaries are located 
at y = ±12. 8A, and the Z boundaries are located at z = ±12. 8A. We take the relativistic 
Harris model as the initial current sheet. The magnetic field, the density of plasmas and 
their distribution functions are described by B — Botanh{z/\)x, n{z) = uq cosh^^(2;/A) and 
f± oc n{z) exp[— r^{£ — P±Uy}/T], respectively. In the above equations, Bq is the magnitude 
of magnetic field in the lobe (background region), uq is the number density of plasmas in 
the current sheet, (3± — v±/c are the drift velocities for each species; /9+ = +/3 for positrons, 
/3_ = — /3 for electrons. is the Lorentz factor for P [F/3 = [1 — e is the particle 

energy, u is the relativistic four velocity of w = [1 — (f /c)^]~^/^ ■ v and T is the temperature. 
In the first simulation, we set T = mc^ and (3 = 0.3. We set no driving force to excite 
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instabilities, so that instabilities arise from thermal noise. We set no background plasmas. 
The total energy is conserved within an error of 0.3% throughout the simulation run. 

3. Results 

Snapshots at three characteristic stages of the simulation are presented in Figure 1. The 
three left panels in Figure 1 show color contours of the plasma density at (a) t/Tc — 46.0, 
(6) 64.0 and (c) 82.0, where time is normahzed by the light transit time Tc = A/c. The three 
middle panels show color contours of the electric field Ey at the corresponding time; the 
white lines show contours of B^- The three panels show the particles' energy spectra for the 
whole simulation. The horizontal axis shows the particles' energy normalized by mc?. 

As shown in the top three panels in Fig. 1, we observe a kink-type modulation of the 
current sheet. We confirmed that this is a signature of the RDKI, a relativistic extension of 
the DKI, which is a long- wavelength, current-driven instability in a thin current sheet. In 
the top middle panel, one can see the typical structure of Ey components of the polarization 
of electric fields of RDKI. The sign of Ey is positive in the yellow regions, while it is negative 
in the blue regions. Two types of regions appear alternatively along the current sheet and 
they are anti-symmetric with the neutral plane {z — 0). Note that the Ex B direction using 
the lobe magnetic fields and the observed sign of Ey is consistent with the Z-displacement 
of the plasma bulk motion. The symmetric E^ components are also observed and their 
values are a little less than or comparable to those of Ey. The E^ components are negligible. 
Since the system is uniform in the X direction, there is no electrostatic E^ fields. After 
t/Tc ~ 50, the instability turns into its non linear stage. The three middle panels in Fig. 1 
are snapshots of the current sheet at t/r^ = 64.0. The current sheet is strongly folded within 
a zone of —4 < z/X < 4. Due to the Z-displacement by RDKI, red (originally yellow) 
regions, where Ey > 0, are placed in a row around the neutral plane {z = 0) in the midmost 
panel of Fig. 1. The total kinetic energy in the system has increased by 30% at this time. 
Importantly, we observe a clear sign of the particle acceleration in a high-energy tail in the 
energy spectrum of the middle-right panel in Fig. 1. We will discuss the acceleration process 
later. The kinked-shape of the current sheet is not stable, because fragments of the anti- 
parallel currents (±Jz) pull each other. As they start to collide, the system turns into the 
"mixed" stage. The snapshots at this stage are presented in the bottom three panels in Fig. 1 
at t/Tc = 82.0. Due to the enhanced diffusion by sheet collisions, the total kinetic energy 
has increased by 170% at this time. In the energy spectrum of the bottom-right panel, one 
can recognize both a remnant of the non-thermal tail and global plasma heating. Finally, 
the mixed current sheet slowly evolves into the broadened current sheet, which is 3-4 times 
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thicker than the initial state. The saturated level of the total kinetic energy is 330%-340% 
at t/Tc — 200 or later. The spectrum looks unchanged from the mixed stage; one can still 
observe the non thermal tail. 



4. Discussion 

Now, we again take a look at the instabihty during its hnear stage. The observed growth 
rate [uJi/Q^c) is plotted as a function of ky\ in Fig. 2, where ky is the wavenumber and 
is the gyro frequency Vtc = ^dl = (eB) / {'-fmc'^) . The growth rate and wave length of the 
most dominant mode of RDKI are Ui/Qc ~ 0.035 and kyX ~ 1.0, respectively. We have also 
calculated eigen functions and their growth rates of the RDKI, by linearizing relativistic two- 
fluid equations. The theoretical growth rates for the present case arc shown as a solid line 
in Fig. 2. The simulation results are consistent with them except for shorter wavelength of 
kyX > 1, where the instabihty is suppressed by kinetic effects. We observe a shght signature 
of the relativistic drift sausage instabihty (RDSI) along with RDKI. However, RDSI is less 
influential on the current sheet deformation and particle acceleration. Note that no velocity- 
shear-driven mode such as KHI is excited, ant that is because there is no velocity shear 
between drifting plasmas in the current sheet and empty background plasmas. Also note 
that LHDI is not excited in an e"^ case, because there is no delay of a positron's response to 
an electron's. 

Next, we discuss an acceleration mechanism that is found at the non linear stage of 
RDKI. By analyzing trajectories of highly accelerated particles, we found that the main 
site of particle acceleration is the central channel around z ^ 0, where the current sheet 
was originally located. For simple discussion, we call this channel the "acceleration channel 
(AC)". We illustrate the acceleration mechanism in schematic views in Fig. 3. In the upper 
view, the structure of the polarization of Ey fields at the hnear stage are illustrated. They 
are anti symmetric with respect to Z and both positive and negative regions are alternatively 
located along the sheet. Passing through their meandering orbits, some particles can resonate 
with the Ey fields. Some of them gain their energy and others lose their energy. At the non- 
linear stage, due to the Z-displacement of the current sheet, the positive Ey regions are 
now located around the AC. nsidc the positive Ey regions, typical values of the fields are 

~ 0.6 — 0.7-Bo and Ey ^ 0.2 — 0.3-Bo- Contrary to the reconnection (Zenitani & Hoshino 
2001; Jaroschck ct al. 2004) in which > 1 in an acceleration site, the value of 

is less than unity inside the positive Ey regions, and so low-energy particles travel into the 
E X B direction. However, fewer populations of high-energy particles can travel across the 
positive Ey regions, because their Larmor radii are larger than the half-wavelength of 
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RDKI (n/ky). They gain their energy from the electric field Ey, cross the current sheet and 
move into the neighboring positive Ey regions. Since becomes larger as particles gain 
more energy (r^ oc e^^"^ in non-relativistic case and oc e in relativistic case), accelerated 
particles usually satisfies the crossing condition of > n/ky in the following positive Ey 
regions. In this way, they are successively accelerated in the AC, passing through multiple 
positive Ey regions into the ±Y direction. The above condition can be interpreted as the 
"trapping condition" inside the AC. Higher accelerated particles tend to stay in the AC, 
because larger prevents them from escaping. This is why the non-thermal tail is enhanced 
in the energy spectrum. We think this type of enhanced acceleration can be found in any 
kink- type instability in a thin current sheet, as long as its wavelength is comparable to a 
particle's Larmor radius. For example, it may be applicable to acceleration of ions, when a 
thin current sheet of ion-electron plasmas is modulated by KHl. The folded structure evolves 
into the mixed stage when its Z-displacement (A2;) is nearly the same as the half-wavelength 
of RDKI {Az n/ky vrA). Assuming that Vz is typically Vexb ~ 0.3c in the positive Ey 
regions, the time scale of the non-linear stage {tn) is approximated by rjv ~ ttA/w^ ~ 5 — IOTc. 

After the system turns into the mixed stage, a few number of high-energy particles 
and a lot of low-energy particles in the current sheet interact with the positive Ey regions. 
Then the magnetic fields are dissipated in a time-scale of tens of Tc- Roughly speaking, 
the magnetic field energy is consumed by Joule heating of the Ey fields and the total cur- 
rent Ty in the broadened current sheet. The magnetic energy stored between the folded 
curves is approximated by {7r/ky)BQ/S7r. We assume that the average value of the elec- 
tric fields Ey ~ 0.5Ey — 0.5 x (0.2 — O.S)Bq, where the factor 0.5 represents the fact that 
positive Ey is not uniform in the AC. The total current ly is represented by the zeroth or- 
der current Ty ~ A Jo = cBq/Atv. Thus the dissipation time scale (td) is approximated by 
td ~ {7r/ky){B^/87r){XJoEy)~' - 10 - 20re. 

Next, we briefiy compare RDKI and magnetic reconnection (Zenitani & Hoshino 2001). 
In Table 1, the parameters for several simulation runs of RDKI and for relevant runs of 
reconnection are presented. Throughout the simulation runs, we calculate the ratio of the 
non-thermal kinetic energy to the total kinetic energy (Knonth/K) and their maximum values 
in time are presented in Table 1. In case of run D3, Knonth is not resolved by our method. 
The maximum energies {Smax) are also presented in Table 1. In the case of RDKI, the 
highest energy particles run in the AC into the ±Y direction with the light speed during 
the non-thermal stage and the mixed stage. Their maximum energies are estimated by 
^est ~ ^maxo + GcEy{TN + T£)) ~ s^axo + ^uJcTcTTic^ and the estimated values agree with the 
simulation results. In three cases of reconnection, the evolution seems to be restricted by the 
periodic system size of L^, because two outflow jets travel into the ±X direction. Due to this 
restriction, the upper limits of particle energy are estimated by Sest — ^BqLx/2 in Table 1 



-6- 



in reconnection cases. However, the obtained maximum energies (emax) by reconnection are 
far larger than the maximum energies by RDKI. Prom the viewpoint of field dissipation (e.g. 
striped pulsar winds), RDKI can dissipate more magnetic energy. The linear growth rates of 
the energy conversion {Tcdln{AK) / dt) arc presented in Tabic 1. Generally speaking, RDKI 
grows faster in a relativistic current sheet (T/mc^ > 1) and it converts 2.5 — 3.5 times as 
much energy as the initial kinetic energy in the current sheet in a time scale of td- Moreover, 
RDKI spontaneously occurs everywhere in the X direction, while single reconnection requires 
a wide spatial range in the X direction in order for the outflow jets to escape. Therefore, 
even in runs R3 and D3, in which RDKI grows slower, reconnection takes twice as long to 
dissipate the same amount of the magnetic energy as RDKI. On the contrary, reconnection 
is the more favored candidate as the origin of non-thermal particles. It produces more 
non-thermal energy of Kno^th than RDKI as a result of the high ratio of Knonth/K. 

Our initial results of a three-dimensional simulation show that RDKI occurs almost 
uniformly in the current sheet. Also, Jaroschck ct al. (2004) have presented their three- 
dimensional results of triggered reconnection. They have referred to the importance of 
RDKI in the three-dimensional case and a signature of a two-dimensional-like RDKI is 
observed outside the reconnection region. In case of multiple current sheets structure, such 
as striped pulsar winds, there are several possibilities of enhanced field dissipation. If the 
distance between two currents is less than the critical distance of 2Ac, the two current sheets 
collide with each other and 90% or more of the magnetic energy will be dissipated into the 
particle energy. In typical cases, Ac ~ 7-8, but its dependency on T or /3 should be further 
investigated. We also note that the secondary RDKI takes place in the broadened current 
sheet in some occasions. 

Let us summarize this Letter. First, we investigate the RDKI in a thin current sheet. 
Second, we find a new acceleration process that uses a folded structure of a current sheet. 
The unique point of this process is that AC electric fields evolve into a DC acceleration 
channel. Third, due to the enhanced acceleration, a non-thermal particles are produced. 
Fourth, comparison with reconnection shows that RDKI is more favorable to dissipate the 
magnetic energy in the current sheet structure. We believe that RDKI provides a crucial 
clue to the energy dissipation problem in the current sheets in the universe, including striped 
pulsar winds. 

This work was supported by the facilitates of JAXA and the Solar- Terrestrial Environ- 
ment Laboratory, Nagoya University. 
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Fig. 1. — Left three panels: Color contours of the plasma density at I/tc =46.0, 64.0 and 
82.0. Middle three panels: Ey at the corresponding time; the white lines shows contours of 
Bx- Right three panels: Energy spectra. 
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(a) Linear Stage energy loss 




Fig. 3. — Schematic view of the acceleration mechanism. Once the current sheet is folded at 
the non linear stage, regions with the positive Ey sign are located in a central acceleration 
channel. 
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Run 


T 


Lx 


Ly 




Knonth/ K 


Tcdln(AK)/dt 


^max 


^est 


Dl 


16 


1.0 


25.6 


51.2 


4.7 X 10-2 


2.3 X 10-1 


7.0 X 102 


7.1 X 102 


D2 


1 


1.0 


25.6 


51.2 


7.2 X 10-2 


2.3 X 10-1 


4.8 X lOi 


4.2 X lOi 


D3 


1/16 


1.0 


25.6 


51.2 




1.9 X 10-2 


4.0 


3.0 


Rl 


16 


102.4 


1.0 


51.2 


> 2.5 X 10-1 


6.6 X 10-2 


> 3.6 X 10^ 


5.6 X 10^ 


R2 


1 


102.4 


1.0 


51.2 


> 2.6 X 10-1 


6.5 X 10-2 


> 1.2 X 102 


3.5 X 102 


R3 


1/16 


102.4 


1.0 


51.2 


> 1.1 X 10-1 


7.1 X 10-2 


> 5.3 


2.2 X lOi 



Tabic 1: Three runs for RDKI (Dl-3) and relevant runs of relativistic reconnection (Rl- 
3) are presented. The temperature {T/mc^), the system size [L^, Ly, Lz in unit of A), the 
maximum non-thermal ratio of the kinetic energy [Knonth/ K), the energy conversion rate 
{Tcdln{AK) / dt) , the maximum energy {smax/^c^) and their estimated values {eest/^c^) are 
presented. 



